Single-mode delay time statistics for scattering by a chaotic cavity 
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Abstract 



We investigate the low-frequency dynamics for transmission or reflection of a 
^ ■ wave by a cavity with chaotic scattering. We compute the probability distri- 

bution of the phase derivative 4>' = d(/)/du; of the scattered wave amplitude, 
known as the single-mode delay time. In the case of a cavity connected to 
two single-mode waveguides we find a marked distinction between detection 
in transmission and in reflection: The distribution P{(j)') vanishes for negative 
(j)' in the flrst case but not in the second case. 
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I. INTRODUCTION 

Microwave cavities have proven to be a good experimental testing ground for theories of 
chaotic scattering |]l|. Much work has been done on static scattering properties, but recently 
dynamic features have been measured as well 0. A key dynamical observable, introduced 
by Genack and coworkers |0-§[, is the frequency derivative 0' = d^/dcu of the phase of the 
wave amplitude measured in a single speckle of the transmitted or reflected wave. Because 
one speckle corresponds to one element of the scattering matrix, and because cj)' has the 
dimension of time, this quantity is called the single-channel or single-mode delay time. It is 
a linear superposition of the Wigner-Smith delay times introduced in nuclear physics ^^■ 

The probability distribution of the Wigner-Smith delay times for scattering by a chaotic 
cavity is known [^]. The purpose of this paper is to derive from that the distribution P(0') 
of the single-mode delay time. The calculation follows closely our previous calculation of 
P{4>') for reflection from a disordered waveguide in the localized regime ||^. The absence of 
localization in a chaotic cavity is a significant simplification. For a small number of modes 
A^ connecting the cavity to the outside we can calculate P{(f)') exactly, while for A^ ^ 1 
we can use perturbation theory in 1/A^. The large-A^ distribution has the same form as 
that following from diffusion theory in a disordered waveguide [^^, but for small A^ the 
distribution is qualitatively different. In particular, there is a marked distinction between 
the distribution in transmission and in reflection. 



II. FORMULATION OF THE PROBLEM 

The geometry studied is shown schematically in Figure |1[ It consists of an A^-mode 
waveguide connected at one end to a chaotic cavity. Reflections at the connection between 
waveguide and cavity are neglected (ideal impedance matching). The A^ modes may be 
divided among different waveguides, for example, N = 2 could refer to two single-mode 
waveguides. The cavity may contain a ferri-magnetic element as in Refs. ||TU|,ITT|, in which 



case time-reversal symmetry is broken. The symmetry index (3 = 1 (2) indicates the presence 
(absence) of time-reversal symmetry. We assume a single polarization for simplicity, as in 
the microwave experiments in a planar cavity 0. 

The dynamical observable is the correlator p of an element of the scattering matrix S{u) 
at two nearby frequencies, 

p = Snm{uJ+lSuj)S*^^{uJ-l6uj). (l) 

The indices n and m indicate the detected outgoing mode and the incident mode, respec- 
tively. The single-mode delay time 0' is defined by ||^-|^ 

0' = lim ^, (2) 

with I = \Snm{^)\'^ the intensity of the scattered wave in mode n for unit incident intensity 
in mode m. If we write the scattering amplitude Snm = v/e*'^ in terms of amplitude and 
phase, then 0' = d(j)/duj. We will investigate the distribution of 0' in an ensemble of chaotic 






FIG. 1. Sketch of a chaotic cavity coupled to N propagating modes via one or more waveguides. 
The shape of the cavity is the quartered Sinai bilhard used in recent microwave experiments [0]- 



cavities having slightly different shape, at a given mean frequency interval A between the 
cavity modes. For notational convenience, we choose units of time such that 2tt/A = 1. 

The single-mode delay times are linearly related to the Wigner-Smith 0,^ delay times 
Ti, r2, ..., Tat, which are the eigenvalues of the matrix 



Q = -iS''-— = t/Miag(ri, ..., tn)U. 
auj 



(3) 



To see this, we first expand the scattering matrix linearly in 6uj, 

S{uj ± l6uj) = V^U ± li6ujV^dia.g{Ti, ..., tn)U. (4) 

Since 5" is symmetric for (3=1, one then has V = U. For (3 = 2, V and U are statistically 
independent. Combination of Eqs. (|I|), (^, and (^ leads to [Q 
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The distribution of the Wigner-Smith delay times for a chaotic cavity was calculated in 
Ref. [P|. It is a Laguerre ensemble in the rates /ij = 1/tj, 



P{iii,...,Hn) oc Jll/^i 
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\P^^k)e{^^u). (7) 

for X < 0. It follows from Eq. (|^ that 



The step function 9{x) = 1 for x > and 9{x) 
{J^i^i) = 1; ci result that was known previously 

To calculate the joint distribution P(/, 0') from Eq. (||), we also need the distribution of 
the coefficients Ui and fj. This follows from the Wigner conjecture |jl3[|, proven in Ref. [Q, 
according to which the matrices U and V are uniformly distributed in the unitary group. 



The calculation for small A^ is now a straightforward integration, see Sec. |T|. For large N 
we can use perturbation theory, see Sec. |V|. 

Because of the uniform distribution of U and V, independent of the r^'s, we can evaluate 
the average of (p' directly for any A^, 
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We define the rescaled variable 0' = 0'/(0') = Ncj)', that we will use in the next sections. 



III. SMALL NUMBER OF MODES 

For N = 1 there is no difference between the Wigner-Smith delay time and the single- 



mode delay time. In that case 1 = 1 and 0' = 0' is distributed according to |T^ . [T5[] 

P(0O = C/30' ^' exp(-i/3/0O^(0O- (9) 

The normalization coefficient C/3 equals (27r)^^/^ for /3 = 1 and 1 for (3 = 2. 

Now we turn to the case N = 2. By writing out the summation in Eq. (^ for / and 0', 
one obtains 0' = r+ + ar_ with t± = |(ri ± T2) and 

/Il2ll2ill2ll2i H= =f=iH= * /ir\\ 

= |-ui| |t;i| + \U2\ \V2\ + U1U2V1V2 + U1U2V1V2, (lOj 

a = (Kph|2-|«2nt;2p)//. (11) 

To find the joint distribution P{I, a) we parametrize U in terms of 4 independent angles, 

jT _ ( cos7exp(— mi) sin7exp(— iai — 'iQ;2) \ /.„x 

y — sin7exp(— iaa + ia2) cos7exp(— zas) y' ^ '^ 

with Oj G (0, 27r) and 7 G (0,7r/2). The invariant measure d/i oc |Det fi'|d7nidaj in the 
unitary group follows from the metric tensor gf, defined by 

Ti dUdU^ = ^gijdxi(lxj, {xi} = {7, ai, ^2, as}- (13) 

hi 

The result is 

dyU oc sin27 d7 JJdaj. (14) 

i 

The joint distribution function P{t^, t_) follows from Eq. (^. For f3 = I one has 

P{r^, r_) = ^|r_|(r^ - r!)-^exp {-r^{rl - r!)-^) ^(r+ - |r_|), (15) 

while for (3 = 2 

Pin,r-) = Irlirl - t'_)-<^ exp (-2r+(r^ - rl)-') ^(r+ - |r_|). (16) 

First we consider the case (3 = 1, n ^ m. Because of the unitarity of U, one has 
l^iP = |^2p and |f2p = |miP- Therefore a = and 0' = r+, so 0' is independent of /. 
Integration of Eq. (^) over r_ results in 

P(0') = f 0'"'(0'' + 20' + 2) exp(-2/0')^(0')- (17) 

In this case (as well as in the case A^ = 1), 0' can take on only positive values, but this is 
atypical, as we will see shortly. From Eqs. (^^) and (|12]) we find the relation between / and 
the parametrization of U, 



/ = sin^ 27sin^(Q;3 — ai — a2)- (18) 

The distribution of I resulting from the measure ([T^) is 

p{i) = \i-y^e{i)e{i-i), (19) 



in agreement with Refs. ||16| , p!7| . 

For the case N = 2, f3 = 1, n = m we use that Ui = Vi, U2 = V2 and obtain the 
parametrization 

1=1 — sin^ 27 sin^(a3 — ai — 0:2), (20) 

a = (cos 27)//. (21) 



The distribution P{I, a) resulting from the measure (|T4D is 



P(/, a) = — /i/2(l - /)-i/'(l - Ia^)-'^^e{I)e{l - I)e{l - la^). (22) 

27r 

The joint distribution of / and (p' = 2(f)' takes the form 
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P(/,0') = / dr_ / dr+ P{T^,T_)P(l,a = ^^^_I±\J__ (33) 

The distribution of / following from integration of P{I, a) over a is given by Eq. ([ToD with 
/ — i> 1 — /, as it should. The integrations over t^,t^, and J, needed to obtain P(0') can be 
evaluated numerically, see Fig. ^. Notice that P{4>') has a tail towards negative values of 0'. 
For A^ = 2,/5 = 2 it doesn't matter whether n and m are equal or not. Parametrization 
of both U and V leads to 

/ = (1 - xi) (1 - X2) + X1X2 + 2y (1 - a;i) (1 - X2)xiX2 cos r], (24) 

a={l-xi-X2)/I, (25) 

with a measure d/i oc da;ida;2d?7 and xi, 0:2 G (0, 1), r/ G (0, vr). The joint distribution P{I, a) 
is now given by 

P{I, a) = \l^/^e{I)e{l - 1)6(1 - la"). (26) 

Integration over a leads to [P^PT[ 

p{i) = e{i)e{i-i). (27) 

The distribution P(/, </>') follows upon insertion of Eqs. ( |16|) and (|26|) into Eq. (|23D . Numer- 
ical integration yields the distribution P(0') plotted in Fig. ^ As in the previous case, there 
is a tail towards negative cf)'. 



1^ 




FIG. 2. Distribution of the single- mode delay time in the case of preserved time-reversal sym- 
metry (/? = 1). The curves for A^ = 1,2 follow from Eqs. (0), (p^, and (p3|). The curve for 
A^ ^ 1 follows from Eq. (|36|), and is the same for n = m and n ^ m. The delay time (j)' = (j)' /{(/)') 
is rescaled such that the mean is 1 for all curves. Data points are the result of a Monte Carlo 
calculation in the Laguerre ensemble (with N = 400, n ^ m representing the large- A^ limit). 
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FIG. 3. Same as in Fig. |2|, for broken time-reversal symmetry (/? = 2). The curves for A^ = 1, 2 
and for A^ ^ 1 follow from Eqs. @, (p^, and (|36|). There is no difference between n = m and 
n ^ m for any A^. The large A^-result for /3 = 2 is the same as for P = 1. 



IV. LARGE NUMBER OF MODES 



We now calculate the joint distribution P{I,(j)') for A^ ^ 1. First the case n ^ m will 
be considered, when there is no distinction between (3=1 and (3 = 2. In the large A^-limit 
the vectors u and v become uncorrelated and their elements become independent Gaussian 
numbers with zero mean and variance 1/A^. We first average over v, following Ref. |p. We 
introduce the weighted delay time W = !((>'. The Fourier transform of P{I, W) is given by 
x{p, q) = {exp{ipl + iqW)). The average over v is a Gaussian integration, that gives 



xip,q) = {det{l-tH/Nr'), 



*,-,Tn 



TT * T I 1 /-* T I *- 

H = pu u + ^qy^ u + u u 

where Mj = UiTi. The matrix H has only two nonzero eigenvalues. 



A± = |( g5i + p ± \/2pqBi + q^B2 + p^ 



Bk 



EklV- 
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Performing the inverse Fourier transforms and returning to the variables 0' and / leads to 



P{I,^') = (Ar3//7r)^/2exp(-iV/) ( {B2 - Bl)-^^^ exp 



-NI- 



B2 — Bi 



9{I). (32) 



The averages over Ui and Tj still have to be performed. 

Up to now the derivation is the same as for the disordered waveguide in the localized 
regime 0, the only difference being the different distribution of the Wigner-Smith delay 
times Tj. The absence of localization in a chaotic cavity greatly simplifies the subsequent 
calculation in our present case. While in the localized waveguide anomalously large r^'s lead 
to large fluctuations in Bi and B2, in the chaotic cavity the term /if ' in Eq. (0) suppresses 
large delay times. Fluctuations in B^ are smaller than the mean by a factor 1/viV. For 
A^ :» 1 we may therefore replace Bk in Eq. ( ^2]) by (Bk). 

To calculate the average of Bi and B2 we need the density p{t) = (X^j S{t — Ti)) of the 
delay times. It is given by [H 



N 



pM = 7r-i^{r+-T){r-T^), 



27TT 
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(33) 



for T inside the interval {t^,Tj^). The density is zero outside this interval. The resulting 
averages are (Bi) = N~^ and (^2) = 2N~'^, which leads to 



P(J,0O = {N^I/TTy/^exp(^~NI 



1 + 



9{I). 



(Recall that 0' = 0'/(0') = N((>'.) Integration over 0' or / gives 

P{I) = Nexp{-NI)e{I), 



P(0O 



1 + 



'3/2 



(34) 

(35) 
(36) 



This distribution of / and 0' has the same form as that of a disordered waveguide in the 
diffusive regime |^|[. 

We next turn to the case n = m and f3 = 1. (For (3 = 2 there is no difference between 
n = m and n ^ m.) Since Ui = Vi in Eq. (|^) we have 

/=|Cor, 0' = Re^, C, = Y.rtul (37) 

The joint distribution P(Co,Ci) has the Fourier transform 

x(Po,Pi, go, gi) = (exp(zpoReCo + igoImCo + ipiReCi + iqilmCi)). (38) 

Averaging over u we find 

x{Po,Pi,qo,qi) = {exp{-L)), (39) 

L = i ^ In [l + iV-2(po + p,T,y + N'\qo + qiT.f] . (40) 

i 

Fluctuations in L are smaller than the average by a factor 1/A^. We may therefore approx- 
imate (exp(— L)) ^ exp(— L). Because N~'^{po + piTiY + N~'^{qo + giTj)^ is of order 1/A^, 



we may expand the logarithm in Eq. (40). The average follows upon integration with the 
density (PBD, 

Inverse Fourier transformation gives 

P{Co,C^) = —-^exp{-N\Co\' - '^N^\C,\' + N^Re CqC*). (42) 

The resulting distribution of 0' and / is 

P(/, 0') = (iV3//27r)i/2 exp (-^A^/ [l + (0' - 1)^] ) 6(1). (43) 



It is the same as the distribution (|3^) for n ^ m, apart from the rescaling of / by a factor 
of 2 as a result of coherent backscattering. 

The distribution (^) of 0' for A^ ^ 1 is included in Figs. ^ and ^ for comparison with 
the small A^-results. 



V. NUMERICAL CHECK 

We can check our analytical calculations by performing a Monte Carlo average over the 
Laguerre ensemble for the r^'s and the unitary group for the Mj's and f j's. For the average over 
the unitary group we generate a large number of complex Hermitian N x N matrices H. The 
real and imaginary part of the off-diagonal elements are independently Gaussian distributed 
with zero mean and unit variance. The real diagonal elements are independently Gaussian 

8 



distributed with zero mean and variance 2. We diagonalize H , order the eigenvalues from 
large to small, and multiply the n-th normalized eigenvector by a random phase factor e*"" , 
with a„ chosen uniformly from (0,27r). The resulting matrix of eigenvectors is uniformly 
distributed in the unitary group. 

The Laguerre ensemble (|^) for the rates ^i = I/tj can be generated by a random matrix 



of the Wishart type |]T8|,|T9[. Consider a A^ x {2N — 1 + 2//5) matrix X, where X is real 
for /? = 1 and complex for /5 = 2. (The matrix X is neither symmetric nor Hermitian.) 
The matrix elements are Gaussian distributed with zero mean and variance (|a;„mP) = 1- 
The joint probability distribution of the eigenvalues of the matrix XX^ is then given by 
Eq. (j^). The results of our numerical check are included in Figs. ^ and 0. The large-X limit 
is represented hj N = 400, n ^ m. The analytical curves agree well with the numerical 
data. 



VI. CONCLUSION 

We have investigated the statistics of the single-mode delay time 0' for chaotic scattering. 
For a large number N of scattering channels the distribution has the same form as for 
diffusive scattering [Q0, but for small N the distribution is different. The case N = 2 is 
of particular interest, representing a cavity connected to two single- mode waveguides. For 
preserved time-reversal symmetry and detection in transmission {(5 = l,n ^ m),we find that 
(f)' can take on only positive values, similarly to the Wigner-Smith delay times. In contrast, 
for detection in reflection (or for broken time-reversal symmetry) the distribution acquires 
a tail towards negative 0'. These theoretical predictions are amenable to experimental test 
in the microwave cavities of current interest [^ . 
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